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We report on the charge carrier dynamics of superconducting titanium nitride (TiN) in the fre- 
quency range 90-510 GHz (3 - 17 cm"^). The experiments were perfomed on a 18 nm thick TiN 
film with a critical temperature of = 3.4 K. Measurements were carried out from room tempera- 
ture down to 2 K, and in magnetic fields up to B = 7 T. We extract the real and imaginary parts 
of the complex conductivity as a function of frequency and temperature, directly providing the 
superconducting energy gap 2A. Further analysis yields the superconducting London penetration 
depth Al- The findings as well as the normal state properties strongly suggest conventional BCS 
superconductivity, underlined by the ratio 2A(0)/A:_b Tc = 3.44. Detailed analysis of the charge 
carrier dynamics of the silicon substrate is also discussed. 



For the last few years promising research has been 
done on highly yet homogenously disordered TiN super- 
conducting thin films, since this material has been iden- 
tified to be a prime candidate for the superconductor- 
insulator transition^ (SIT). Besides fundamental physi- 
cal interest, TiN has recently gained attention as a ma- 
terial for microwave resonators^. Today, one can access 
a broad range of experimental data displaying the dc- 
properties of TiN superconducting thin films, e.g. com- 
prehensive dc-transport experiments or local scanning 
tunneling spectroscopy at millikelvin temperatures and 
high magnetic fields. Very intriguing phenomena like a 
energy gap*^ for T > Tc which might serve as key for un- 
derstanding the pseudogap^ in high- Tc superconductors, 
voltage threshold behaviour^ or a notable peak in the 
magnetoresistance-'- have been observed at the verge of 
the SIT. All these phenomena, as well as the very nature 
of the SIT itself, are not yet entirely understood. For the 
last years, several theoretical approaches were studied®. 
One common belief is that superconductivity is destroyed 
on a macroscopic level when entering the normal state. 
However, microscopic islands embedded into an insulat- 
ing or metallic matrix remain superconducting but have 
lost any mutual phase correlation. This picture is sup- 
ported experimentally by magnetoresistance effects^, lo- 
cal STM-, voltage threshold behaviour^'^ and numeri- 
cally by Monte-Carlo simulations There is, however, 
very little known about the electrodynamic response of 
superconducting thin films, which are known to show the 
SIT. Infrared (IR) reflectivity measurements^^ were per- 
formed on several TiN thin film samples with different 
degrees of disorder and broadwave microwave corbino 
spectroscopy on various InO thin film samplesi^. With 
our work, we extend the studied frequency range well 
below the far IR down to 3 cm^^. This range provides 



direct access to the superconducting properties, e.g. the 
real and imaginary parts of the complex conductivity, the 
energy gap and the London penetration depth as func- 
tions of frequency and temperature below and above the 
energy gap. Identical experiments were performed with 
two similar TiN films. Since both TiN films lead to the 
same results, we will just discuss one of them. The pre- 
sented data was taken on a 5 x 5 mm^ large TiN film 
of 18 nm thickness grown via plasma enhanced atomic 
layer deposition at T = 400 °C on a 0.73 mm thick sil- 
icon (110) substrate covered with a 10 nm thick Si02 
film. The critical temperature is = 3.4 K and the 
sheet resistance i?n(300K) = 91.6 Vl. Considering these 
values, we expect a TiN film, which is at the boundary 
between bulk and quasi-two-dimensional superconductiv- 
ity. Fig. [T}d displays the sheet resistance i?n of the TiN 
film as a function of temperature. The data was taken 
in a standard 4-point-measurement using Van-der-Pauw 
analysis^"^. Starting at room temperature, R\2 decreases 
linearly from 91.6 VL to 85 Q.. Immediately before the 
superconducting transition, a small peak shows up at 15 
K, see the inset in Fig. lb. This peak emerges from the 
competition between weak localization, electron-electron 
interaction and supercondcuting fluctuations. In TiN 
fllms with a higher degree of disorder this peak is much 
more pronounced-- while it is absent in clean films. Since 
the absolute value of quantum contributions to conduc- 
tivity is much smaller than conductivity itself in what 
follows we will not take it into account. At 4 K, i?g de- 
creases rapidly and drops to immeasurably small values 
at Tc = 3.4 K. 

We performed transmission and phaseshift measure- 
ments covering the spectral range 3-17 cm~^ for 2-292 
K and magnetic fields up to 7 T. Tuneable backward- 
wave-oscillators (BWO) emitted coherent radiation of 
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FIG. 1. (color online) (a) Real part of the substrate conduc- 
tivity versus temperature for different frequencies from 4-23 
cm~^. Below 50 K, the conductivity drops to zero and the 
substrate remains transparent (b) Sheet resistance of the TiN 
film versus temperature. We emphasize the good agreement 
between dc and THz approach (solid line and open circles re- 
spectively). The Inset shows a close up of the sheet resistance. 
The small peak is a signature of quantum contributions. 
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FIG. 2. (color online) (a) Transmission Tf and (b) phaseshift 
over frequency (/it/tj versus frequency at 2 K and 10 K below 
and above Tc = 3.4 K. The solid line is a calculation according 
to BCS theory. The change in both properties indicate a 
strong frequency dependent complex conductivity. The inset 
shows the temperature evolution of the reduced energy gap 
together with the BCS prediction leading to 2A(0) = 8.14 
cm~^. 



high power, which was then detected by a Golay cell. The 
TiN film sample and a bare Si substrate sample as refer- 
ence were mounted into an optical ^He bath cryostat ori- 
ented as shown in the inset of Fig. 2a. Both transmission 
and phaseshift were measured with the same quasiopti- 
cal Mach-Zehnder interferometer—"—. Transmission Tf 
and phaseshift over frequency (pt/tJ were fitted simulta- 
neously using well known equations based on Fresnels 
formulas for multiple reflections^^. Fitting parameters 
were the real and imaginary parts of the refractive in- 
dex, n and k, respectively. Considering these two quan- 
tities, one obtains any optical function desired, e.g. the 
complex conductivity a as follows 



y/n + ik = ei(a;)+ie2(w) = 1 (T2{uj)-\ cri(w) (1) 

The open circles in Fig. [TJd is the sheet resistance^^ 
i?n = (dtTo)~^ (d = 18 nm fllm thickness) extracted from 
the THz data. We obtain the same results for Rfj within 
an accuracy of a few percent for both dc and THz ap- 
proach. Fig. la shows real parts of the complex con- 
ductivity of the silicon substrate versus temperature for 
several frequencies between 4-23 cm~^. Below 50 K, ai 
drops to zero and remains there down to our base tem- 
perature of 2 K. This implies that the substrate is purely 
dielectric in the superconducting regime of TiN. A de- 
tailed discussion of the THz properties of the substrate 
is provided in the appendix. Fig. 2 displays Tf and 
(/>f/w versus frequency for 2 K in the superconducting 
and 10 K in the normal state. The solid line is a cal- 



culation according to the BCS theory applied to both 
Tf and ^t/w simultaneously, with n and k expressed by 
CTi and (72 of the Mattis-Bardeen formalism22i2i (dirty 
limit). The oscillations are caused by interference due to 
multiple reflections inside the substrate, which behaves 
like a nonsymmetric Fabry-Perot resonator. At 10 K, the 
Fabry-Perot oscillations show a uniform pattern, which 
does not vary with frequency. As implied by the sheet 
resistance in Fig. lb, this is due to the dispersionless 
metallic conductivity of the TiN film at frequencies much 
lower than the scattering rate. As soon as the critical 
temperature is passed, Tf and (/it/a; change dramatically, 
and Tf rises well above the normal conducting values at 
2 K and for frequencies larger than 6 cm~^. This can be 
understood in the two-fluid model; the number of quasi 
particles contributing to ai is reduced due to the forma- 
tion of Cooper pairs, which do not contribute to (Ti for 
flnite frequencies below 2A. For frequencies lower than 6 
cm~^, Tf drops below the normal state Tf because of the 
increasing reflection since the superconducting film acts 
like a dielectric mirror, (pt/io, however, decreases for all 
frequencies. The considerable changes below Tc allude a 
strong frequency dependent complex conductivity. 

From the raw data we extracted the superconduct- 
ing properties of the TiN film: we perform single-peak 
fits for each Fabry-Perot oscillation separately. For each 
peak, i.e. each resonance frequency, we obtain a pair 
of ei and €2, which is then used to calculate the com- 
plex conductivity. Real and imaginary parts of (7(0;) are 
shown in Fig. 3 for several temperatures below and above 
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FIG. 4. (color online) Energy gap versus magnetic field (open 
circles) together with the prediction based on the phenomeno- 
logical two-fiuid model. The fit leads to _Bc = 5 T. 



FIG. 3. (color online) (a) Real and (b) imaginary parts of 
the complex conductivity versus frequency for four different 
temperatures below Tc = 3.4 K and one well above. The 
solid lines are fits according to BCS theory applied to both 
(Ti and (T2 simultaneously. The kink in a\ indicates the energy 
gap, which shifts towards lower frequencies upon increasing 
temperature. The insets show the London penetration depth 
versus frequency and temperature, respectively. The solid line 
in the upper inset is the two-fiuid model prediction yielding 
Al(0) = 730 nm. 



Tc. The solid lines are fits according to Mattis-Bardeen 
theorji2£i2i (dirty limit) applied to both o\ and 02 si- 
multaneously. As can be seen, u\ and 02 are well de- 
scribed by the BCS theory. Starting from 2 K, a\ is 
strongly suppressed for frequencies around 7 cm~^ and 
approaches the normal state frequency independent con- 
ductivity (To = 6350±400 il^^cm^^. The minimum in a\ 
indicates to the energy gap that gets less pronounced and 
shifted towards lower frequencies upon increasing tem- 
perature. This trend is displayed in the inset of Fig. 2b 
together with the BCS expressions^ 



^=tanh(^^| 
A(0) 1 T A(0) / 



(2) 



The fit leads to a zero temperature gap of 2A(0) — 8.14 
cm"-'^ and consequently a BCS ratio of 2A(0)/fcsrc = 
3.44 very close to the weak-coupling prediction^ of 3.53. 
Upon cooling, G2 rises notably from zero and approaches 
(72 oc w"^, which is consistent with a\ cx 5(a; = 0) at 
zero temperature according to Kramers-Kronig relations. 
This trend is also seen in the frequency dependence of the 
London penetration depth Al, which was calculated from 
(72 viapi^ 



Al = 



47rW(72 



Assuming the above mentioned a; "^-behavior, Ai gets 
frequency independent. In contrast to T < Tc, one finds 
a nearly constant A^ at 2 K, which underlines the theoret- 
ical prediction. For T < Tc, however, Al shows a strong 
rise towards high frequencies. We attribute these high 
penetration depths to the increasing number of quasi- 
particles due to thermal and photon induced Cooper- 
pair breaking. Upon increasing photon energy hw > 2 A, 
Cooper pairs are destroyed, and hence the superconduct- 
ing state is weakend. Extrapolating the Ai(w) curves to 
zero frequency leads to the data shown in the inset of Fig. 
3a, where the squared ratio Al(0)/Al(T) versus temper- 
ature is shown. The data was fitted within the two-fluid 
model^^ 



Ai(0) ^._(T_ 



(4) 



The fit yields the London penetration depth Al(0) = 
730 ± 50 nm at zero temperature. 

The data presented so far strongly supports the inter- 
pretation as a conventional BCS superconductor. We 
also performed experiments in magnetic fields up to 
_B = 7 T at 2 K utilizing a superconducting cryomag- 
net optical access and studied the influence on the en- 
ergy gap. The magnetic field was applied parallel to the 
film. The energy gap was extracted by fitting Tf and 
(jytl^^ simultaneously for each magnetic field. The fitting 
parameters were 2A and Tc while the other material pa- 
rameters, plasma frequency ojp and scattering rate F, 
utilzed for the fit were calculated separately from the 
normal state at 4 K and were taken as constants. Fig. 4 
displays energy gap versus magnetic field together with 
a fit via the phenomenological approach^— 



(3) 



A(i?) _ L_(B_ 
A(0) V Uc 

yielding a (upper) critical field of Be 



(5) 

5 T. The gap 
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decreases monotonically with increasing magnetic fiek 
until the pair-breaking hmit is reached and superconduc 
tivity is totaUy suppressed. The analysis of the raw date 
for higher fields did not indicate an energy gap. Th( 
zero-field value of 2A of the different setups are in gooc 
agreement (7.5 cm^^ and 7.2 cm^^). 

In summary, we have studied the charge carrier dy 
namics of superconducting TiN films with a low degret 
of disorder at temperatures well below = 3.4 K anc 
frequencies above and below 2A. The superconducting 
properties are well described within the BCS theory. Fu- 
ture efforts will have to extend the experimentally acces 
sible range towards lower temperatures to shed light oi 
the yet not understood phenomena close to the SIT. 
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FIG. 5. (color online) Ts of silicon substrate versus temper- 
ature and frequency. The frequency dependence shows the 
well developed Fabry-Perot oscillations, that extend over the 
entire temperature and spectral range except for the low-Ts 
region at low frequencies around 100 K. 



Appendix: THz properties of the Si substrate 

Since our measurements of TiN thin films were carried 
out as combined measurements of film and substrate, we 
had to analyze the properties of the bare substrate in 
order to disentagle TiN film and Si substrate contribu- 
tions. We therefore measured a 1 x Icm^ slab (thickness 
0.78 mm) of the same type of doped silicon, onto which 
our films are deposited. We found our raw data to be 
in excellent agreement with previous THz-time-domain 
measurements^^ of doped silicon, and thus used the same 
way of analysis. In Fig. 5, versus temperature and fre- 
quency is displayed. Starting at room temperature, Ts is 
about 0.4, slightly decreasing towards lower frequencies. 
Upon cooling, Tg is reduced and decreases to zero at low 
frequencies at around 100 K. Ts rises again upon further 
cooling, and reaches unity at 50 K for each frequency. 
We do not observe changes of the Fabry-Perot oscilla- 
tions between 50-2 K, which indicates a purely dielectric 
Si substrate. 

Fig. 6c displays cri versus temperature and frequency. 
For a given T > 50 K, cti decreases monotonically to- 
wards higher frequencies. The slope is maximal at around 
100 K. From 100-50 K, the dispersion dies off rapidly and 
we do not observe any residual conductivity for lower 
temperatures. Our single-peak fitting routine does not 
assume any particular electronic model, i.e. the extracted 
complex conductivity can be utilzed to identify the most 
suited theory describing the charge carrier dynamics of 
our substrate. Fig. 7 shows cti and (T2 versus frequency for 
three different temperatures together with fits according 
to the Drude modeli^ 
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FIG. 6. (color online) Charge carrier dynamics extracted from 
Ts. (a) Charge carrier concentration (solid line: Fermi-Dirac 
distribution) (b) mobility and scattering rate are determined 
by fitting o\ within the Drude model, (c) real part of con- 
ductivity (Ti versus frequency and temperature. The data is 
extracted from Ts by single-peak fitting. 



cri(a;, T) -(- ia2{uj, T) = 



.ao(T)c^F(T)- 



l-cj2F(r)-2 1 



a;2r(T)-2 
(A.l) 

where cto — NTme ^ is the dc conductivityi^ (e elemen- 
tary charge, m electron mass) and F = is the scat- 
tering rate (inverse scattering time), both used as fitting 
parameters. Clearly, the Drude model is well suited to 
describe the charge carrier dynamics, however, small de- 
viations of tT2 are apparent. We assume that the devia- 
tions are caused by Ts being close to zero and therefore 
not showing pronounced Fabry-Perot oscillations. This 
does not affect cti, since this quantity is related to the 
amplitude of T^. However, it affects (T2, since this quan- 
tity depends on the position of the oscillations with re- 
spect to the frequency, which can not be assigned with 
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FIG. 7. (color online) Real part a\ and imaginary part 02 
(inset) of the complex conductivity versus frequency for three 
temperatures. Solid lines are fits according to eqs. A.l & A. 2. 



sufficient accuracy here. Utilizing the Drude fits of cti, 
we calculated the charge carrier concentration N shown 
in Fig. 6a. The carrier concentration N is supposed to 



satisfy the Fermi-Dirac distribution that reads=i 

with Nd the donor density, Ny=6 the number of equiv- 
alent valleys for Si, m„ = 0.32rno the effective mass, h 
Planck constant, Boltzmann constant and AEiy the 
donor ionization energy. The solid line in Fig. 6a is 
the Fermi-Dirac distribution with AEn = 62 meV and 
Nd — 3.5 X 10^''^ cm~'^. Localization of the charge carriers 
explains the vanishing conductivity at low temperatures 
as soon as a certain activation temperature is gone be- 
low. Fig. 6b shows the mobility ^ and the scattering 
rate-- F = e(m/i)~^ extracted from the Drude fits. With 
increasing temperature, the charge carrier movement is 
constrained more and more due to the enhanced phonon 
scattering, leading to a reduced mobility and increasing 
scattering rate. 

In conclusion, our Si substrate shows charge carrier 
dynamics, which can be assigned to the doping. Above 
an activation temperature of about 50 K, the frequency 
dependence of the complex conductivity is well fitted by 
the Drude model. Below 50 K, the Si substrate is purely 
dielectric, which enables us to unambigously disentangle 
Si substrate and TiN film contributions. 
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